Fractal Behavior of the Shortest Path Between Two Lines in Percolation Systems 
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Using Monte-Carlo simulations, we determine the scaling form for the probability distribution 
of the shortest path, I, between two lines in a 3-dimensional percolation system at criticality; the 
two hues can have arbitrary positions, orientations and lengths. We find that the probability dis- 
tributions can exhibit up to four distinct power law regimes (separated by cross-over regimes) with 
exponents depending on the relative orientations of the lines. We explain this rich fractal behavior 
with scaling arguments. 



I. INTRODUCTION 

There has been considerable recent activity ana- 
lyzing P{t\r), the probability distribution for the length 
of shortest path, between two points separated by Eu- 
clidean distance r in a percolation system H-^] . This pa- 
per extends that work by determining the scaling form 
of the distribution of shortest paths between two lines 
of arbitrary position, relative orientation and lengths in 
3-dimensional systems. These configurations are impor- 
tant because they much more accurately model the con- 
figurations used in oil recovery in which fluid is injected 
in one well (one of the lines in our configuration) and 
oil is recovered at a second well (the second line in our 
configuration); the wells may, in reality, be at arbitrary 
orientation and of different lengths. 

The scaling form for the 2-points configuration 
in which the 2 points arc located at {{L — 
r)/2, L/2, L/2), {{L + r)/2, L/2, L/2) in a system of side 
L has been found to be |§ 
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The exponents 



and ip are universal and the 



constants a and b depend on lattice type. In 3D, the val- 
ues of these exponents have previously been found to be 
§ 5£ = 2.3 ± 0.1, dmin = 1.39 ± 0.05, = 2.1 ± 0.5, and 
i/j ~ 2.5 ± 0.5. The first stretched exponential function, 
/i , reflects the fact that the shortest path must always be 
at least equal to the distance r between the two points; 
the second stretched exponential function, /2 , reflects the 
fact that the lengths of the shortest paths are bounded 
because of the finite size, L, of the system. 



We find that the scaling form for the 2-lines configu- 
ration has the same form as that found for the shortest 
path distribution between two points with the exceptions 
that: (i) the power law regime of the distribution as rep- 
resented by the term (.£/r'^min)-s« in Eq. (|^) is replaced 
with up to four different power law regimes (separated 
by cross-over regimes) with exponents depending on the 
relative orientations of the lines and (ii) the Euclidean 
distance, r, in Eq. (^ between the two points is replaced 
by the shortest Euclidean distance between the two lines. 
The lengths of the lines affect the sizes of the power law 
regimes. 



II. SIMULATIONS 

Monte Carlo simulations were performed using the 
Leath method and growing clusters from 2 sets of seeds — 
one for each line. The length of the shortest path between 
the two lines is the sum of the chemical distances from 
each set of seed sites to the point where a cluster started 
at one set of seeds meets a cluster started from the other 
set of seeds Q. The cluster growth for a given real- 
ization is terminated when the two clusters meet. For 
parallel line configurations, in which the probability dis- 
tributions decay rapidly, we use the method of Ref. 0] 
to obtain good statistics for shortest paths that have very 
low probabilities. We use the memory management tech- 
nique described in pT[ | to perform simulations in which 
the growing clusters never hit a boundary of the system. 

The clusters that are created and included in our anal- 
ysis are of all sizes, not just the incipient infinite cluster. 



III. NON-PARALLEL WELLS 
A. Co-Planar Lines 

1. Equal Length Symmetric Lines 

We start by considering relatively simple configura- 
tions of the type shown in Fig. |l](a) in which the lines 
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are co-planar, of equal length and are positioned sym- 
metrically. We will study configurations in which the 
lines are of unequal length [see Fig. 0(b)] and/or are not 
positioned symmetrically [see Fig. jjc)] in the following 
sections. In all of these configurations, r is the shortest 
Euclidean distance between the two lines. 

Figure || contains log-log plots of P{l\r)^ the shortest 
path distribution for ?' = 8 and various values of 9. We 
have chosen r = 8 so that the initial cutoff is present; 
for smaller values of r, lattice effects destroy this initial 
cutoff. Since the focus of this paper is the power law 
regimes, not the initial or final cutoffs, in all later figures 
we will choose configurations with r = 1 so that the ex- 
tent of the power law regimes is as long as possible. The 
exception to this will be cases in which 6 is very small 
where small r introduces other lattice effects. 

We note that after the initial peak in each distribution, 
there is a power-law regime, the slope of which, gi{0), in- 
creases with increasing 6. We will call this power law 
regime the "2-lines regime." Simple scaling arguments 
imply that if the lengths of the lines were infinite, these 
2-lines regimes would continue indefinitely. For finite line 
lengths, we would expect that, for large £, the distribu- 
tions would exhibit a crossover to a power law regime 
with the same exponent as that for a configuration with 
two points — because for large £ the long paths travel far 
enough away from the lines that they appear to be points. 
In this regime, the power law exponent has the value of 
that of two points, 2.35 For the plots in Fig. ||, in 
order to see the power regimes clearly, we have chosen 
the lengths of the lines long enough that this crossover 
occurs after the maximum value of £ in the plots. 

In Fig. ^, we plot ge{0) versus 9. The plot suggests 
that ge{9) diverges as 9 goes to zero. We attempt to fit 
this function with a power of l/sin{(9)/2) and find the 
best fit with the function 



sm{e/2) 



0.4 • 



(4) 



This form is not based on any fundamental theory; it sim- 
ply has the properties that f{9 = 180°) = gt{9 = 180°), 
it diverges as 9 goes to zero and it fits the intermediate 
points reasonably well. 

The crossover between the 2-line regime and the 2- 
point regime is illustrated in Fig. ^ where in each panel 
we plot P{£\r) for fixed 9, and various values of W. As 
expected, the larger the length of the lines, the higher the 
value of £, the value of £ at which the crossover occurs. 
Quantitative analysis of the crossover behavior is given 
in Section III.C. 



2. Point-Line Configurations 

We next study configurations in which one line has 
zero length (i.e., a point) and the other is a line of fi- 
nite length. This is the extreme case of the configura- 
tion in which the two lines have different length. We 



will study the case where both lines have finite length 
in the next section. We in fact study the three configu- 
rations shown in Figs. |^(a)-(c). The plots of P{£\r) for 
these configurations are shown in Fig. The plots have 
a power- law regime with exponent —1.75 for the config- 
urations of Fig. H(a) and Fig. ||(c), and exponent —2.2 
for the configuration of Fig mh). We denote this regime 
the "point-line regime." FigT]^ shows the crossover from 
point-line behavior to 2-points behavior. 



3. Unequal Length Symmetric Lines 

We can now study configurations of the type shown in 
Fig. 1(b) in which the lines are of different lengths, Wi 
and W2 ■ For such a configuration we would expect three 
power-law regimes: (i) for small £ such that the two lines 
appear to be infinite, a 2-line regime, with slope depen- 
dent on 9, (ii) a point-line regime, with slope —1.75, for 
values of £ large enough that the shorter line appears to 
be a point, and (iii) a 2-point regime, for even larger val- 
ues of £ where both lines appear to be points. Plots of 
P{£\r) for such configurations arc shown in Fig. ^ and 
are consistent with our expectations. 



4- Complex Configurations (Unequal Length Non-Symmetric 

Lines) 

The last of the co-planar configurations is of the type 
shown in Fig. 0(c). In general, based on the reasoning 
above, for configurations of this type we would expect 
P{£\r) to have four power-law regimes. For the configu- 
ration shown in Fig. 0(c), in which Wia ^ W2 <€. Wu, 
the power law regimes would be as follows: (i) a power 
law regime corresponding to the angle 9 between the seg- 
ments Wia and W2, (ii) a power law regime corresponding 
to the angle tt — 9 between segments W2 and Wib, (iii) 
a point-line power law regime entered when £ 3> Wib, 
and (iv) the 2-points regime. Figure || is a plot of P{£\r) 
which shows this behavior. 



B. Non-Coplanar Lines 

For non-coplanar lines, for I 3> r, the fact that the 
lines are not co-planar should be irrelevant; what is rel- 
evant is the effective angle between the lines. This angle 
is obtained by sliding the lines toward each other along 
the line of shortest Euclidean distance between the lines 
(without changing their orientations) until they touch; 
the lines are then coplanar and the angle between them 
determines the behavior of P{£\r). Figure |l^ contains 
plots for two configurations which illustrate this: (i) two 
coplanar lines with r — 1, 9 = 90°, and W = 256, and 
(ii) the same configuration with the second line trans- 
lated out of the plane by distance 8. We see that while 
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there is some difference in the plots for small /, the slope 
of the 2-lines regime is the same for the two plots. 



C. Scaling of the Crossover Between Power Law 
Regimes 

We define the value of ^ = ^ at which P{i\r) crosses 
over from one power-law regime to another power-law 
regime as the value of £ where straight lines fit to the 
power law regimes A)etween which the crossover takes 
place, cross. In Eq.(P the values of £ at which the lower 
and upper cutoffs occur scale independently as r'^min 
and L'^min^ respectively. By extension, we would ex- 
pect that all characteristic values of the distribution, in- 
cluding crossovers between different power-law regimes, 
would also scale as X'^min where X is the length in the 
system which controls the crossover. Thus, in analogy 
with the scaling of the most probable value of £, £* , 



I* 



(5) 



we would, in fact, expect that the value of £, £ at which 
P{£\r) crosses over from 2-lines behavior to 2-points be- 
havior scales as 



where 



r + 2M^sin(6l/2) 



(6) 



(7) 



is the maximum Euclidean distance between the two lines 
and ci{0) is a slowly varying function of 9. In order for a 
2-lines regime to exist, the 2-lines regime cutoff^ must be 
greater than £*, the maximum value of the distribution. 
That is, 



ci[r -I- 2M^sin(6l/2)]'^min > cr'^mm, 



which implies 



W > 



(c/ci) 



2sin(6'/2) 



(9) 



In Figs, ^(a), (b), and (c), the insets contain plots of £ 
vs. Tmax- For 9 = S° , the scaling exponent is consistent 
with Eq.(^ but for 9 — 29° and 180 the scaling exponent 
is 1.0 ±0.1. 

Using the same reasoning which led to Eq. (|), we 
would expect the crossover from point-line to 2-points 
behavior to scale as 



(10) 



because W is the length which controls this crossover; 
as seen in Fig. 0, the larger the value of W, the larger 
the value of £ at which the crossover from point-line to 
2-points behavior occurs. However, as seen in the inset 



in Fig. 0, the crossover length scales with an exponent 
1.0 ±0.1 not d,„i„. 

Finally, we would expect that for different length lines, 
the crossover from 2-lines behavior to 2-point behavior 
would scale as 



£ = C3(9)W^ 



{W2 < Wi] 



(11) 



because W2 is the length which controls this crossover; 
as seen in Fig. ^, the larger the value of W2, the larger 
the value of £ at which the crossover occurs. Again, the 
insets in Fig. || indicate that the crossover scales with 
exponent 1.0 ± 0.1. 

We cannot explain why sometimes the crossover scales 
with dmin and sometimes it scales with an exponent about 
1. It is, of course, possible that corrections-to-scaling are 
strong and that we are not seeing the true asymptotic 
behavior of the scaling of the crossover. If this is the 
case, the question still remains as to why the corrections 
to scaling are strong in some configurations and not in 
others. This area is a subject for further study. 



IV. PARALLEL WELLS 



A. Simple Configurations 

As with non-parallel wells we first consider the simple 
configurations shown in Fig. |ll|(a) in which the parallel 
wells are of the same length. Figure |2|(a) plots P{£\r) vs 

for r = 1 and various W. The initial decay of the plots 
increases with increasing W because the longer the wells, 
the lower the probability for long shortest paths. Eventu- 
ally, all plots cross over to a power-law regime with slope 
consistent with that for two points. To see if this initial 
decay is a lattice effect. Fig. O plots of the scaled distri- 
butions r'^minP(^£/r'^min\W) for various r and W where 



(8) the aspect ratio. 



R = 



W 



(12) 



is fixed at i? = 32. Changing r and W but keeping R 
fixed results in scaling all lengths in the geometry by the 
same factor and the plots collapse as expected. Again, 
we note the initial strong decay of the distribution fol- 
lowed by a 2-point power-law regime. The good collapse 
for small x = l/r''"'" indicates that the strong initial 
decay is not a lattice effect. 

Because of the small values of £ at which the crossover 
to the 2-point regime occurs it is difficult to differentiate 
between a power-law and (stretched) exponential decay. 
We will proceed as if the decay were either a power law 
with slope g{R) or equivalently an exponential with "ef- 
fective slope" g{R)- 

One might argue as follows that the initial decay for 
power parallel lines must be exponential: Since the 2- 
lines regime of the probability distribution for a parallel 
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well configuration must always decay faster than the 2- 
lines regime of a configuration with small but non-zero 
and since we believe gi{9) goes to infinity as 9 goes 
to zero, the decay for parallel lines must be exponential 
(i.e., faster than any power law decay). This, however, 
need not be the case. In order for a 2-lines regime to 
exist, Eq. (||) must hold. So as we decrease 9, we must 
increase W , increasing the aspect ratio R, to maintain a 
2-lines regime. But since the effective slope for parallel 
wells, g{R) increases with increasing i?, the decay can be 
a power law and still always have a greater slope than 
the configuration with small but non-zero 9. 



B. Complex Configurations 

The treatment of the more complex configurations 
shown in Figs. |l^(b) and (c) follows that of non-parallel 
wells. P{tj for configurations of the type in Fig. pT|(b) 
would contain an initial 2-line regime with slope g{R — 
W2/r), a point-line regime, and finally a two-point 
regime. P{i) for configurations of the type in Fig. [ll|(c), 
with Wib ^ W2 ■€i Wia would contain an initial 2- 
line regime with slope g{R — Wib/r), a 2-line regime 
with slope gi{9 = tt), a point-line regime, and a 2-point 
regime. 



C. Quasi-Euclidean Regime 

When the length of the wells is very large and the dis- 
tance between the wells is small the behavior of the most 
probable shortest path between the wells is the same as 
in a Euclidean space where p=l and all bonds are occu- 
pied. This can be seen in Fig. |lj in which we plot i*, 
the most probable value of the shortest path, versus r 
for various lengths W. For long enough wells, there is a 
regime of r in which 



The shortest path will exhibit Euclidean behavior, i.e., 
£* = r when P{r, W) is of the order unity. Setting 
P{r*,W) = a in Eq. (||), we find 



(13) 



as one would expect in Euclidean space in which the 
shortest path is a straight line path of occupied bonds. 
As also seen in Fig for a given well length, as r in- 
creases, there is a value of r, r* , at which the behavior 
crosses over to that of 3D percolation. We can develop a 
simple expression for r* as follows: The probability that 
all bonds in a straight line path between two wells sep- 
arated by distance r are occupied is p^. The probability 
that one or more bonds in the straight line path is not oc- 
cupied is thus l—Pc and the probability that one or more 
bonds in the W straight line paths between the wells are 
not occupied is (1 — p^)^ . The probability that at least 
one straight line path has all bonds occupied is then 



ln(l - gi/^) 
Inpc 



(15) 



P(r, W) = 1-{1~p: 



(14) 



In Fig. |T5| we plot the observed values of r* and the val- 
ues predicted by Eq. ([l5| ) with a value of a = 0.55 which 
gives the best fit to the observed values. 



V. RELATIONSHIP BETWEEN PARALLEL 
WELLS AND "CLOSE TO PARALLEL" WELLS 

For a given r, we expect that a configuration with small 
but non-zero angle will have a power-law regime slope 
very close to the (effective) power-law regime slope of a 
configuration of parallel lines with the same W. This at 
first leads to a seeming paradox: if we increase or decrease 
W, but keep the angle of the non-parallel wells fixed, the 
slope of the 2-lines regime of the non-parallel well config- 
uration doesn't change as discussed in Section III.A.l. 
However, if we consider this configuration as a paral- 
lel configuration, changing W changes the aspect ratio 
which changes the power-law regime slope as discussed 
in Section IV. A. This seeming inconsistency is resolved 
as follows: on the one hand, for a 2-lines regime to exist, 
W must be at least as large as the value given by Eq. (||. 
If W is too small, there will be no 2-lines regime and 
both the parallel and small angle configurations will look 
like the configuration for 2 points. On the other hand, 
if W is increased, keeping r and 9 fixed, the greater the 
deviation from parallel lines and there is no reason why 
the parallel and small-angle configurations should have 
the same slopes in their power-law regimes. 

Thus, only for the very small range of W for which 
the power-law regime exists and for which the configura- 
tion with small but non-zero 9 is "close to parallel" (i.e. 
the difference between the values of r and r^ax is small) 
should the slopes of the parallel configuration and the 
configuration with small but non-zero 9 be equal. That 
is, 

g[W/r] « ge{9), (16) 
where g{R) is defined in Section IV. A. 

VI. DISCUSSION AND SUMMARY 

Motivated by the need to more realistically model the 
geometries found in oil recovery activities, we have de- 
termined the scaling form for the distribution of shortest 
paths between two lines in 3 dimensional percolation sys- 
tems. Using simple scaling arguments we explained the 
rich fractal behavior of the shortest path in these sys- 
tems. A number of open questions, however, remain: 
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(i) From first principles, can one develop an expres- 
sion for gi{9)l An exact expression for gi for two 
points in 2-dimensions was obtained by ZifF § us- 
ing conformal invariance arguments. Possibly this 
approach could be extended to find gi for point-line 
and 2-line configurations, at least in 2-dimensions. 

(ii) How is the fact that the crossover from one power- 
law regime to another does not scale with the ex- 
ponent dmin explained? Is this simply an artifact of 
corrections-to-scaling which would disappear if we 
could simulate much larger systems or is the scal- 
ing of the crossover actually anomalous in certain 
configurations? 
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FIG. 1. Example configurations of two non-parallel lines which are studied, (a) Simple configuration of lines of equal length, 
(b) Configuration of lines of unequal length (Wi > W2). (c) Configuration in which shortest distance between lines does not 
terminate at the ends of lines {Wai < W2 < Wu,). 




10'' 10^ 10^ 

FIG. 2. P{p(r) vs i for configuration of two lines of equal length with r = d>, = (from bottom to top) 3° (filled square), 
6°, 12°, 20°, 40°, and 180° (unfilled square). The corresponding well lengths W are 4890, 2445, 1224, 737, 374 and 128, 
respectively. The plots are normalized such that the initial sections of plots are coincident. 
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10^ 10^ 

FIG. 4. P{l\r) vs (, for configuration of two lines of equal length, (a) r = 8, 6 = 3°, = (from top to bottom) 38, 76, i52, 
304, 1216, and 2432 {h) r = \, 9 = 29°, W = (from bottom to top) 8, 17, 33, 66, 132, 265, and 529, (c) r = 1, = 180°, 
W = (from bottom to top) 8, 16, 32, 64, and 128. For all plots, the larger the value of W, the larger the value of £ at which 
behavior changes from 2-lines behavior to 2-points behavior for which the slope is —2.35. The insets plot the crossover value, 

£, vs. Tmax. 
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(a) 



(b) 



(c) 



W 



w 



r 

• * 



FIG. 5. Example configurations in which one line is of finite length W and one is of zero length (i.e., a point). In all cases, 
the shortest distance between the point and the line is r. 
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FIG. 6. P{(.\r) vs £ for configuration of a point and a line with r — 1 and W = 128. From top to bottom, the plots are 
for the configurations shown in Figs.^(a), (c), and (b) respectively. We see that the slopes in configurations where the point 
is closest to the end of the line [Fig. g (a) and (c)] are the same (with some initial difference) and they are different from the 
slope in the configurations in which the point is closest to the middle of the line (Fig. Hb). 
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10^ 10^ 

FIG. 7. P{i\r) vs i for configuration of a point and a line in which the point is closest to the end of the line with r = 1 and 
W = (from bottom to top) 2, 4, 8, 16, 32, 64, and 128. For all plots, the larger the value of W , the larger the value of £ = £ at 
which the behavior changes from point-line behavior to 2-points behavior. The inset plots I vs. W . 
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10'' 10^ 10^ 

FIG. 8. P(£[r) vs I for configurations of two lines of different lengtlis witli r = 1 and Wi = 128. (a) 6 = 7°, W2 = (from top 
to bottom) 16, 32, and 64. Three power law regimes can be seen: the 2-lines regime, the point-line regime and the 2 points 
regime, (b) 6 = 180°, W2 = (from bottom to top) 4, 8, 16, 32, 64, and 128. Only the first two power law regimes can be seen: 
the 2-lines regime and the point-line regime (the 2-points regime would require even larger values of t) . 
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10^ 10^ 10"* 

FIG. 9. P(£\r) vs I for configurations of two lines of different lengtlis which "overlap" [see Fig. 0(c)] with 9 = 7° , r — 1, 
W\a ~ 32, Wib = 128, and 11^2 ~ 256. Four power law regimes are present: the 2-line (9 — 7°) regime (slope —3.0), the 2- line 
(9 — 180° — 7°) regime (slope —1.2), the point-line regime (slope —1.75), and the 2-points regime (slope —2.35). 
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I I I I ■ ■ ■ ■ I I I I I ■ ■ ■ ■ I I I I I ■ ■ ■ ■ I I I I l_ 

10^ 10^ 10^ 

FIG. 10. P{l\r) vs t for configurations of two fines of equal length. The co-planar configuration has r = 1, 9 = 90°, and 
W = 256 and the lines are co-planar. The non-coplanar configuration is obtained from the co-planar configuration by moving 
one of the lines a distance 8 perpendicular to the plane defined by the coplanar lines. One sees that for large i, the power law 
regimes of the two plots have the same exponent. 
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FIG. 11. Example configurations for parallel wells, (a) Simple configuration of wells of equal length, (b) Configuration of 
wells of unequal length {W\ > W2)- (c) Configuration in which shortest line between end of one well does not terminate at end 
of other well {Wai <W2 < Wu)- 
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\=f Jf^ min 

FIG. 13. P{£\r) vs £ for configurations of two parallel lines of equal length with {W,r) — (from top to bottom) (32,1), (64,2), 
(96,3), (128,4), (160,5), and (196,6) (b) plots of (a) scaled with the variable x = (Jr'^rnm, The plots in (b) collapse nicely as 
would be expected since they all have the same aspect ratio, W/r. The good collapse for small x indicates that the small x 
behavior is not a lattice effect. 
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FIG. 14. Most probable £ vs r for configurations of two parallel lines of equal length, (a) W = 16, r = 1, 2, 4, 8, and 14. 
(b) W = 32, r = 1, 2, 3, 4, 8, and 16. (c) = 64, r = 2, 4, 5, 6, 7, 8, 16, and 32. (d) Combined plot of (a), (b) and (c). The 
upper and lower dashed lines have slope dmin (1.374) and 1.0, respectively. The larger the value of W, the larger the value of 
r at which scaling crosses over from Euclidean behavior to fractal behavior. 
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FIG. 15. Value of r at which behavior changes from Euclidean to fractal, r* , W . The dashed line is a prediction of Eq. 
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